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We apply our recently developed model of a Bose condensate of quantum kink wave in solid 
He4 to understand recent torsional oscillator experimental results of the citical velocities and the 
effect of the steady and oscillating rotations at around 0.1 degree K. When the D.C. rotation is 
present we find a decrease of the Q factor given by Q _1 oc f s f x Q,d/ujto where f s f is the superfluid 
fraction; Qjy, the D. C. angular rotation velocity, u>to, the torsional oscillator oscillating frequency. 
We estimate the AC critical velocity Q A rtt as that required to generate a kink wave of wavevector 
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27T/ 'Ld where L^ is the distance between nodes of the dislocation network. We generalize this to 
include a steady rotation and find a D. C. critical velocity fig' 4 * oc (O^P*) ' 2 . Estimates for both the 
steady and the oscillating critical velocities are in order of magnitude agreement with experimental 
results. We have also examined an alternative mechanism of kink tunnelling through a node in the 
dislocation networm and find that there is also a dependence on the torsional oscillator frequency: 
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^d'* = [^a^^to^] 1 ■ The DC critical velocity Sl c £ lt is ten times higher than the experimental 
Pr value. 



PACS numbers: 67.80.-s 



Since the discovery of an increase (1 per cent) in the solid 4 He moment of inertia in 
torsional oscillator (TO) experiments at around 200_mKJl|, there have been renewed in- 



terests in its low temperature physical properties 



g. 



Many novel physical behavior are 



manifested, such as a very small direct flow and a very small critical velocity (f^ rf « 10 
rad/s). Recently TO experiments are carried out in the presence of both a steady (DC) and 
an oscillating (AC) rotation[8|, |9_|. An increase in damping is observed which increases with 
the DC rotation speed. When the AC rotation velocity is below the critical value, there 
is also a DC critical velocity Q c £^ which is three orders of magnitude larger than the AC 
critical velocity. The DC rotation does not affect the shear modulus. The AC speed of the 
Kubota group [9], 60/im, is higher than Q c J lt . No critical DC velocity was observed. 

We have recently studied the physics of kink waves of dislocations of density n^ and their 
Bose-Einstein condensation (BEC) [10]. The BEC of the kinks makes possible dissipationless 
movement of the dislocation lines. The motion of a dislocation corresponds in part to a 
circular motion of many He4 atoms, each by a different amount. An estimate of the fraction 
of He4 atoms can be obtained by weighting with respect to the strains. With this the 
corresponding " superfluid fraction" due to the motion of the dislocation lines is found to be of 
the order of nd(iQL m , a magnitude that is consistent with current experimental results. Here 
ao is the lattice constant, L m is the mosiac size. The dislocation motion does not produce any 
net linear motion of the atoms and thus will not generate any direct superflow. In this paper 
we apply our model to understand the experimental results of the citical velocities and the 
effect of the DC and AC rotations. We estimate the AC critical velocity VL c ^ lt as that required 
to generate a kink wave of wavevector 2ti /Ld where Ld is the distance between nodes of the 
dislocation network. When the DC rotation is present we find a decrease of the Q factor 
given by Q^ 1 oc f s f x Q d /uto where f a f is the superfluid fraction; Qd, the D. C. angular 
rotation velocity, uto, the torsional oscillator oscillating frequency. We have examined two 
mechanisms for a DC critical velocity Q c £ lt : (1) The DC rotation generates kinks with a 
time independent displacement. Oscillating kinks are in turn generated from this state by 
the oscillating rotation. (2) Similar to the Josephson effect, the combination of the DC and 
AC rotation can cause a steady current of kinks across nodes of the dislocation network 
when the DC rotation is fast enough. We find that for both mechanisms tt c £ lt oc {Vt c ^ lt ) 1 / 2 
where Q c J lt is the critical AC angular velocity. For the second mechanism, there is also 
a dependence on the torsional oscillator frequency: tt c £ lt = [fi5P*u>ro 27I"] 1 ' 2 . Using current 



experimental estimates for the different physical parameters, we find Q c j[ lt of the same order 
of magnitude as the experimental value. An estimate of the DC critical velocity Vt c ^ lt with 
the first mechanism is of the same order of magnitude as the experimental results; with the 
second mechanism the critical velocity is ten times higher than the experimental value. We 
hope this paper will stimulate further experiments and provide tests of the validity of our 
picture. We now describe our results in detail. 

As is well known 111], for a body of density p rotating with angular frequency Q, the 
quantity of interest is 

F = E-tt-M (1) 

where E, M are the energy and the angular momentum measured with respect to a co- 
ordinate system fixed in space. For example, for a simple rotation at frequency Q, con- 
sider a body rotating at velocity v. Then E[y) = 0.5 J drpv 2 r 2 and M = ujdrpr 2 . 
Minimizing F with respect to v we get v = Q, as we expected. At this frequency 
F{u = tt) = F = -0.5fi 2 f drpr 2 

In the presence of only a time dependenct oscillating rotation with an angular frequency 
^ro{t) — ^A^^-Pi^Tot) caused by a torsional oscillator, if the dislocations also move with 
the entire solid, the energy of the system would be E = 0.5 / drp[flTo(t) x (r + u)] 2 . Here u 
is the displacement due to the dislocations. In our picture, as the He4 is rotated , the kinks 
remain in the zero momentum condensate relative to a space fixed coordinate system. This 
reduces the kinetic energy of the system. We obtain E' = 0.5 J drp[QTo{t) x (r + u + Au)] 2 . 
Here flro{t) x Au is the reduction in the angular velocity where Au comes from the 
motion of the dislocations relative to the rotating solid. The "superfluid fraction" is given 
by f af = (E - E')/E . 

We first estimate the critical velocity when only the time dependent oscillation with 
angular velocity ftTo(t) is present. We consider that a critical velocity is reached when 
it becomes possible to excite a kink wave so that it is possible to lower F; AF becomes 
negative. For a network of dislocations the lowest wavevector is of the order of k$ = 2tt / ' L& 
where Ld is the distance between nodes. This wave vector can be further increased if the 
dislocation moves close to defects (He3) which provide further pinning. In that case, the 
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ocity will become higher. We think a lot of the recently observed hysteretic 



121 ] is related to this issue. 



An example of a kink wave oscillating with frequency ujto is given by: 

\ijj >= sinuTot/2 | — k > + cosuJTot/2 \ko > . (2) 

The average velocity of this state is given by 

v =< i/j\v\iIj >= Vq COSUTot (3) 

where 

v = hko/m*, (4) 

m* is the effective mass of the kink wave. We next proceed to estimate AF. 

If there are N kinks per unit length, each with velocity v the velocity at which the 
dislocation moves out can be estimated as follows. For the dislocation to move out by a 
lattice constant a, each kink has to move a distance of 1/N. The time it takes to do this is 
At = l/(Nv). Denoting the position of a dislocation by c, the speed of the dislocation due 
to the finite speed of the kink is 

dc/dt = aexp(icoTot)/At = aexp(iu:Tot)Nv. (5) 

A dislocation with the Burger's vector along the x direction at some point (a,., c y ) along the 



z axis causes an atom at position (x, y) to move by u x (r — c), u y (r — c) 14j. The atomic 
displacement depends on the position of the dislocation. 

Ou/dt = -V c u • dc/dt. (6) 

is the corresponding velocity of an atom a distance r' away from the moving dislocation due 
to the finite speed of the kinks. In the following, we shall assume that the dislocation moves 
in the direction along the Burger's vector which we take to be the x axis. In general, the axis 
of rotation is not parallel to the axis of the dislocation. With respect to the rotation axis, 
the actual displacement should be Ru(R _1 r) where R is the rotaion matrix that can be 
specified by the Euler angles. We shall assume that that this is the case and for simplicity 
of notation, not displayed this dependence at every step. 

When the kink wave is created, the kinetc energy cost for a segment of the dislocation 
between nodes is given by[13j 

AE = NL d (h 2 kl/2m*) (7) 



where <5u is the displacement caused by this state. The term J drp(du/dt) 2 has already been 



included in the kinetic energy of the kink 15| and thus need not be counted twice. 



The change in the angular momentum due to a change of the state (r >> u) of the kinks 
is given by 

AM « / d 2 rp r x <9u(r')/<9t 

There is another term J d 2 rpr x (J^To(t) x Su) which provides a zero time average to 
Q,To{t)AM and thus will be ignored from now on. From eqs. © and fl6]) AM is of the order 

AM « / d 2 rp r x amNvL d d x ,u(r') (8) 

In general u is a sum of contributions from different dislocations located at different positions 

Cj : u = Xli u o( r ' — c i)- AM can be written as a sum of contributions from each of the 

dislocations. 

AM ~ JZ / d 2 rp(r - cj x amNvL d d x ,u(r' - d) (9) 

t J 

The range of integration of each of these terms is of the order of the mosiac size. Since 
d x iu(v' — Cj) is of the order of l/\r' — q| for a single dislocation, we obatin 

AM « 2ixm*NvL d L 2 Ja (10) 

where L m is the mosiac size. From eq. (7), (1) and the condition that AF = 0, we get the 
critical angular velocity 

Q C X U = AE/AM w ha/{27tm*L d L 2 J. (11) 

Using experimental estimates of L d = 5pm, L m = 20pm, and our estimate m* ~ 0.1mH e 4- 
We get fij*' m 10 _3 /s, of the same order of magnitude as the experimental results. 

We next consider the case where the solid is rotating with a constant angular velocity 
Qd an d ask if it is energetically favorable to start moving the kinks to a state of finite 
momentum. Instead of an "oscillating state" as in eq. (2), we consider the possibility of 
creating simple states \ ± k > . The velocity of the kinks will then just be Vq instead of 
v. Going through the same algebra, we arrive at a DC critical velocity that is the same 
order of magnitude as the AC critical velocity. The experimental DC angular velocity is 
higher than the AC angular velocity by two orders of magnitude. We thus assume a state 
so that the dislocations move with the entire solid with the constant angular velocity ftp. 



In the additional presence of an oscillating driving term so that the total angular velocity is 
Q = Qd + &>To{t), we now consider if the dislocations will exhibit oscillating movements. 
If the kinks do not exhibit the oscillating motion, the kinetc energy saved is given by 

AE = / drpv ■ 5v. 

Here vo = [Hu + ^To(t)] x r is the velocity of the solid in constant rotation. 5v(t) is 
the change in velocity due to the kinks not moving with an oscillating velocity so that 
the core position Sc(t) exhibits a oscillating time dependence relative to the rotating solid. 
The change in velocity now has an additional contribution from the coupling of the steady 
rotation: 

Sv= -V c u- (n D x Ac + dAc/dt). (12) 

dAc/dt w — droit) x c, Ac « — droit) x c/iuro- Because vo is a sum of two terms, 
AE = AEd + AEa contains two contrbutions: those from copupling to Qu and those from 
coupling to Qroit)- The coupling term to the constant DC rotation is given by 

AE D = / drpdr, x r ■ 5v(t) 

which has a zero time average. Thus the DC rotation cannot directly drive the dislocations 
to a finite oscillating velocity. The coupling term to the oscillating rotation is given by 



AE A = / drpd TO (t) x r ■ 8v(t), 



Because Sv is a sum of two terms (eq. [T2|) . AEa = AEai + AEa2 where AEai = 
J drpd TO (t) x r(V c u • dAc/dt), AE A 2 = J drpd TO (t) x r(V c u ■ £l D x Ac). AE A has a 
nonzero time average. The ratio AEai/ E ~ fsf provides for the effective reduction of the 
moment of inertia and is of the order of the "superfluid fraction" f s f. Now Ac and dAc/dt 
in eq. f lT2|) and hence AEai and AEa2 are ninty degree out of phase in time. We thus 
expect AEa2 to provide for a damping term, as is observed in the experiments. The ratio 
Ac/idAc/dt) is of the order of 1/uro, the inverse torsional oscillator vibration frequency. 
The total rotation energy Er of the system is a sum of the rotation energy of the container 
and that of solid He4, Eq. The Q factor is defined with respect to Er- We write E = aEr 
for a constant a. We thus expect the AEa2 term to provide a damping that is of the order 
of Er/Q where 

1/Q w af sf n D /uro- (13) 



Taking a superfluid fraction of the order of 1 per cent, a Qd of the order of 1 rad/s and 
Uto = 2tcx 10 3 rad/s, we obatin an estimate of Q that is of the order of 10 6 a;. Experimentally, 
Q^ 1 ranges from 10~ 6 to 10~ 9 . Our estimate is consistent with this. 1/Q scales with ftp and 
f s f, also consistent with experimental findings. We next examine the critical DC rotation 
field. We have considered two possible mechanisms. We describe them sequentially next. 

(i) We again examine the energetics of creating a kink wave of wavevector 27T /L d . Before 
the kink wave is created, the atoms are at positions V{ + Ui + Au{. Because Aui << Ui, we 
shall neglect the contribution due to Aui below. The angular momentum is now given by 

M = m Y^[r\ + u; + 5ui(t)} x ([ft D + ft TO (t)} x [r ; + u ; + 5ui(t)] + dUi/dt]) 

i 

The velocity Vj is a sum of that due to motion of the kink, duu/dt, and that due to the 
rotation Qp + Qto The corresponding energy is 

E = 0.5m J2 (l fi D + Oto(*)] x [ri + u ; + 5 Ui (t)]) 2 + N L d h 2 / (2mL 2 d ) . 

i 

Recall that before the kink wave is created, the energy is 

E = 0.5m J2 ([^d + Oto(*)] x ( r i + u :]) 2 • 

i 

The change in energy is thus 

AE w NL d h 2 /(2m kink L 2 d ) + m ^(fi D + fi TO ) 2 [0.5^(t) 2 + <&«*(<) fa + «<)]. 

i 

We now look at AF, the change in F as a kink wave is created. 

In general, r >> u(r), after discarding contributions with zero time averages, we obtain 

AF w iVL^V(2m fcmfc L2)-m^[^^(t)r i + 2fi T o(t)fi D r i (5^(t) + fi T o(t)r i ^/9t] (14) 

j 

The last term is the same as in the AC case. Since \Qd\ » \^a\, the term 
2QTo{t)ft>Di r i5uii(t) is much smaller than Q 2 D Suu(t)ri and will be ignored. In this sum there 
is now a new driving term — m ^ Q 2 D 5uu(t)ri that couples to a constant change of position 
of the kinks. Consider, for example, the wave function <f>(z) oc [1 + sm(2-n-z/L d )} which is a 
linear combination of the state \k — > and the states \k = ±27r/L d >. This state has a 
constant shift in the kink position. Once this state is created, the oscillating Hamiltonian 
can couple the states to an oscillating state such as (f)'(z) oc cosojtoI^ + sm(27iz/L d )]. 



The displacement 8uu is of the order Ld, the new term is of the order of magnitude 
—mVL 2 D LdL 2 n / a 2 . Substituting this into eq. f|T4|) and setting AF = we thus arrive at a 
critical DC angular velocity of the order of magnitude 

n% u « [n c fv /L d } 1 ' 2 (15) 

From this we obtain an estimate of f2^ rf of the order of rad/s, the same order of magnitude 
as the expeimental results. 

(ii) We have considered an alternative mechanism due to the onset of the tunnelling of 
a kink wave across the node in the dislocation network. We find a critical angular velocity 
given by 

n crit = [c^Vo27r] 1/2 . (16) 

This critical velocity is of the order of 10 rad/s, a little higher than the experimental value. 
For this mechanism, Vt c ^ lt is a function of the torsional oscillator frequency whereas this is 
not true with the other mechanism. We explain this next. 

We have investigated this by modelling our calculation along the lines similar to the 
Josephson effect with the node of the network modelled as the insulating barrier . Un- 
der the oscillating rotation, due to the centrifugal force there is an effective "poten- 
tial" qAV ~ mL^LdSl 2 /a driving the kinks of the dislocations across the node. Q and 
hence qAV contains both a DC contribution qAVo ~ mL 2 n Ld0£,/a and an AC part 



U_Dc 



qAVA ~ mL^L^Q D Q a cos(u At) /a. As we learned from the Josephson equations |17|], a 
current of kinks can develop across the node that contains a term given by 

J = qAVAsmuTotcos(5o + qAVDt/h)/(htdTo), 

whereto is a constant phase difference. The critical velocity is reached when a DC component 
of the current is developed across the junction. This happens when the quantum energy 
associated with the oscillation frequency ItjjJto is equal to the effective potential applied due 
to the centrifugal force qAVo- We obtain a critical DC angular frequency given by eq. ffT6|) . 
We close this paper with other issues that we have considered. 

As is mentioned above, in general, the axis of rotation is not parallel to the axis of the 
dislocation. The crystal orientation can be specified by two Euler angles (6, $) with respect 
to the rotation axis. (The third angle corresponds to the angle of rotation). The actual 
displacement from the dislocation motion which contributes to the kinetic energy of the 
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particles should be Ru(R _1 r) where R is the rotaion matrix that can be specified by the 
Euler angles. We have explicitly computed this quantity and verified that our results are 
as expected. More precisely we find that J d 2 rr x R<9 x 'u(R _1 r) = 0.5cos2<l>.F(#), F(8) = 
jd 2 r[(l - 2s)(x A cos6-y A cos 3 0) + (cos 2 6 - 1)(3 - 2s) cos 9y 2 x 2 ]/[(s - l)(cos 2 6y 2 + x 2 ) 2 } 
Similarly we obtain J drv ■ R<9 x 'u(R _1 r) = sin(2$)G(#) where G — —[x 2 cos 2 (6) — y 2 ]/[y 2 + 
cos 2 (6)x 2 } 

We were also concerned about possible changes in the phonon dispersion due to the rota- 
tion and its effect on the energetics of the system. We find that the dominant contribution 
to the energy change is given by E 2s d = 0.25(h/N) J2 k . |f2 x e- i \ 2 {2nkj + l)/wfc. where k, j 
specifies the wave vector and branch index of the phonons with frequency Uk, polarization 
e.j and occupation number n^ r Since the phonon frequencies are of the order of 10 12 / sec 
and n is less than rad/s, these corrections are small. 

In summary we apply our recently developed model of a Bose condensate of quantum 
kink wave in solid He4 to understand recent experimental results of the citical velocities 
and the effect of the steady and oscillating rotations. Estimates of the critical velocities and 
the change in the Q value of the trosional oscillator with no adjustable parameters are of 
the same order of magnitude as the experimental results. Their functional dependence on 
system parameters is discussed. We thank Norbert Mulders for helpful discussions. 
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